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Abstract 


A method is presented for the evaluation of the heating efficiency 
required for ignition as a spreading fire closes with fuel. 


An array of thermocouples was implanted in the fuel ahead of the 
fire to obtain the heat absorbed by the fuel prior to ignition. The frac- 
tion absorbed compared to the total that would be absorbed if uniformly 
heated is the effective heating number. The effective heating number 
is represented in its final form as decreasing exponentially with the 
reciprocal of the surface area-to-volume ratio. 


The spreading fire is carried in a regular fuel crib array. Two 
fuel sizes are used in cross-sectional dimension: 0.6 cm. and 1.3 cm. 
Isotherms are generated from the temperature field at ignition and 
presented for both fuel sizes to determine heat absorbed. 


Introduction 


A spreading fire is a series of particle-to-particle ignitions where the energy 
necessary for propagation has its source within the combustion zone where particles 
are undergoing a measurable weight loss through pyrolysis (Fons 1946). These 
consecutive ignitions do not depend on uniformly heated fuel. Frandsen (1971) 
accounted for nonuniform heating ahead of the fire by defining the effective heating 
number. 


Fuel particles ahead of a spreading fire (fig. 1) receive heat from the combustion 
zone, which increases their internal energy. The distribution of heat within similarly 
shaped fuel particles at ignition is dependent on the following factors: location 
of the heat source and heat flow through the surface; surface area-to-volume ratio 
and thermal diffusivity of the fuel. A moisture gradient can also play a major role. 
Assuming that the net heat flux and diffusivity are fairly constant and that the 
moisture gradient is zero, then the surface area-to-volume ratio becomes the signifi- 
cant independent variable. Only in fine fuels (diameter <l mm.) is heating uniform. 
Larger fuels do not achieve uniform internal temperatures because heat impinging on 
the particle cannot diffuse from the surface at a rate comparable to the influx of 
heat. The result is a thermal gradient within the particle that may leave the tempera- 
ture of the central portion, or of one side, relatively unchanged, depending on 
whether the flow of heat is uniform over the surface or unidirectional. 


Although fire spread is attributed to particle-to-particle ignition, the interval 
between ignitions must be averaged over a long time to obtain a rate of spread 
representative of the bulk array. The average rate of spread can then be related to 
the average fuel bed properties. If a constant rate of spread is achieved, the quasi- 
Steady state exists. The fuel bed can then be represented by a minimum unit fuel cell 
that has all the characteristics of the fuel array and retains its usefulness as an 
elemental unit volume in applying the conservation of energy (Frandsen 1971). 
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To employ the concept of the unit fuel cell, the individual particle view must be 
extended to the array within the unit cell. The unit cell then replaces the particle 
as the heat absorber. Consequently, attention is focused on the heat absorbed per 
unit volume, Ppig> rather than on the heat absorbed per unit mass by the individual 


particle where Ph is the bulk density and Qs is the heat of preignition. 


In the unit volume concept the distribution of heat within the particle is less 
important; so the assumption that all heat is uniformly absorbed to a given depth in 
the fuel particle (Thomas 1967) is an equally useful concept. Frandsen (1971) 
extended the individual particle concept to the fuel array by introducing effective 
bulk density. The problem is now resolved to the efficiency of heating the unit fuel 
cell to ignition. The heat per unit volume, Ppig> is necessary to bring the unit 


cell uniformly to ignition. Therefore, the heat per unit volume for ignition at non- 
uniform heating is EPH g» where ce, the effective heating number, is the efficiency. 


The unit fuel cell is not restricted to regular arrays. In field conditions where 
irregular arrays exist, the unit fuel cell must be described by statistical methods. 


The following brief exposition serves to show how e is related to the fire spread 
model developed by Rothermel (1972). 


The rate of fire spread through a fuel array depends on the static fuel parameters, 
Be Pp? and Qi: and the dynamic source function, Ips (Rothermel 1972). The static 


parameters are related to the absorption of energy in the preignition phase, whereas 
the dynamic source function is related to the combustion zone. Their ratio leads to 
the rate of fire spread through the array: 


EI I 
pens ees 


where —€ is the efficiency of converting total reaction intensity, I,, to propagating 
intensity, eS portion of intensity driving the fire (Rothermel 1972). 


€ = Effective heating number 

B = ey on = packing ratio = volume occupied/total bulk volume 

oe Ovendry bulk density 

Pp = Ovendry particle density 

Qe = Heat of preignitaon, the heat. per unit ovendry mass necessary to 


Venite the fuel. 


In this study, the effective heating number was related to the size (surface area-to- 
volume ratio) of the particles making up a regular fuel array of similarly shaped 
particles. “All other variables were held constant. For convenience, the subscript p 
has been dropped from p for the remainder of the paper. Henceforth, the particle 
density will be denoted simply as op. 


Method of Analysis 


Fuel cribs were used in this investigation because of the relative ease of con- 
struction and instrumentation, and their symmetry relative to the fire front. The 
thermal gradient within the particle members of the cell must be obtained to determine 
the heat absorbed up to ignition by the unit cell. 


The unit cell reproduces the crib by consecutively repeating the cell in all 
three cartesian coordinates. Each horizontal layer of cells was alternately shifted 
with respect to adjacent layers to obtain a more uniform array (fig. 2). 


Ignition takes place when a significant portion of the cell undergoes visual 
flame attachment. At ignition, isotherms are constructed from the observed gradient. 
Through symmetry, these isotherms are assumed to form curvilinear surfaces that are 
parallel to the front of the spreading fire. 


The heat absorbed between adjacent isotherms is calculated and summed over the 
area of the isothermal profile to determine the total heat absorbed. The total heat 
absorbed by the unit cell up to ignition is 2V,, where Dis the nonuniform heat absorbed 
per unit volume and V_ is the total volume of the unit cell including the unoccupied 
space. The equivalence of this product to the heat absorbed by the members of the 
unit, celts: 


2V, = LAQ, + tBoQ,. (2) 


Figure 2.--A fuel ertb showing the baste untt cell. At left is the thermocouple arrange- 
ment tn the longttudinal member and above ts the arrangement tn the transverse 
member. The cross-sectional area of a member is A = t2 and the stde area ts 
B= tk, where t ts the member thickness and & ts the spacing on centers between 
fuel members. 


Two identical members of length, 2%, make up the unit cell. They overlap and join at 
right angles (fig. 2). Q is the nonuniform heat absorbed per unit ovendry mass. The 
A and B subscripts refer to the transverse and longitudinal members, respectively. A 
comparison standard for determining the effective heating number was established by 
assuming uniform heating throughout the two members of the unit cell. 


Thomas and Simms (1963) used the moist fuel specific heat and the temperature rise 
to ignition to determine the heat necessary to bring the fuel uniformly to ignition. 
I have employed a method similar to that used by Byram and coworkers (1952), with two 


simplifications: the heat of desorption is ignored! and the water in the fuel is 
assumed to be completely vaporized at 100° C. A temperature-dependent dry fuel specific 
heat (Dunlap 1912), 

Cafe Ce eC, (3) 


is used in place of the constant value used by Byram and coworkers (1952), where: 


C. = 0.266 


2 
iT} 


0.00116. 


Both members of the cell would absorb the same heat per unit mass if the unit cell 
were uniformly heated. For example, the heat absorbed within the transverse member 
would be: 


Hy = Vet£C (100 - T) + fQ + [C, + C, oe + A - Tass (4) 
where: 

Va = Volume of transverse member of unit cell 

f = Fractional fuel moisture content relative to its ovendry mass 

Cc. a Snewiie heat of water 

Cy = Specific heat of ovendry fuel 

ae = Ignition temperature 

7 = Ambient temperature 

Q = Heat of vaporization of water. 


Since Hy = pV Q, and Qn = Q. = Qig for uniform heating, the heat of preignition 


per unit ovendry mass of moist fuel rising from ambient to ignition temperature is: 


G_ = fC, (100 - T,) + #9, + [C 


6 C, oe + aan - ve (5) 

Qj, 1s made up of three products: the heat necessary to bring fuel moisture to 
boiling temperature, the heat of vaporization of fuel moisture, and the heat necessary 
to ignite ovendry fuel. 


The amount of heat nonuniformly absorbed by the unit cell is calculated stepwise 
throughout the fuel member by using isotherms developed from the thermal gradient. 
There are two equations for the heat absorbed by each member within the volume bounded 
by consecutive isotherms: one for T, < 100° C. and one for Bee 100° C., where T,as 
the minimum isotherm. Each equation consists of two parts: 


1. Heat absorbed by all area elements up to T,; 


lHeat of desorption should be included; however, the added refinement is not 
justified when compared to the overall accuracy of these measurements. 


2. The summation of heat absorbed by each elemental area ce or Py bounded by 


ple andedya Upp stor dh. =o.) 
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For example, the heat of preignition for the transverse member is: 


i = 100: 
, Jem Es = 
Cee Qiet one, Halsey ae es oe + Le = he) 
yim 
FS) A. (46.00 = 12), + £0 
A jem+1 J W fe) W 
NC c, (7, + Tz Le =a) (6a) 
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Q, = FIC, (100 - cp) + QO + [C, + Cc, (Tr + T)/2)(T, - i) 
and, 
F : Fie 
j iow j/ 
To< 100) Cie for jp <n 


Ee 100 CY for 5 


W 
r= 


TL. = r100" 4G, tfor jf > mi 


T. = maximum average temperature of Ay for j =n. 


The value for Qp is obtained by applying equations (6a) and (6b) to the longitu- 
dinal member after replacing the summing area A with B. Both results are readily 
evaluated by using a simple iterative computer program to sum the heat absorbed by 
each elemental area. 


Finally, the effective heating number, ce, is evaluated by comparing the ratio of 
the total heat absorbed within the unit volume, ptt (On + Qp), to the total heat that 
would have been absorbed had the unit volume been heated uniformly to ignition, 
2e2t*Q;,. Thus, 


€ = (Qa + Qp)/2Qig- (7) 


Technique 


A unit cell from the upper layer was instrumented to obtain the temperature profile 
within the cell members. Two fuel sizes were used (0.6 cm. (1/4 inch) and 1.3 cm. (1/2 
inch)) in square cross section. Both fuel arrays had a packing ratio of 0.08. The 
packing ratio, 8, is the ratio of the occupied volume to the total volume of the unit 
cell. From the geometry of the bed, 

B = t/2. 5 (C8) 


The horizontal layers lie one upon the other and are alternately shifted, so that 
identical horizontal members are spaced 4t on centers below each other. The greatest 
separation between fuel crib members is in the horizontal direction, the direction of 
fire spread. From equation (8) the spacing, % dimension, is 8 cm. for t = 0.6 cm. and 
16 cm. for t = 1.3 cm. 


Equation (8) is developed from the unit cell shown in figure 2 and applies ex- 
clusively to that arrangement. Intuitively, the rate of spread is dependent on both 
t and 2, but a simple dependency cannot be derived by substituting 8p or (t/2%)p into 
equation (1) for Py> since I, is also dependent on 8 in a less direct manner. 


The instrumented cribs were burned under controlled (no-wind) conditions at the 
Northern Forest Fire Laboratory (Rothermel and Anderson 1966). A uniform fire front 
across the full width of the bed was obtained by igniting alcohol-soaked excelsior in 
the crib spaces at one end. 


Temperature fields within the members of the unit cell were obtained from 5 mil. 
chromel-alumel thermocouples within the longitudinal and transverse fuel members 
(figs. 3 and 4). (The view shown above the fuel array in figure 2 illustrates how 
the transverse member was sectioned to implant the array of thermocouples. Resorcinol 
glue was used to hold the parts together.) The leads immediately adjacent to the 
junction of the thermocouple are parallel to the axis of the transverse member for at 


Figure 3.--Instrumented transverse and longitudinal fuel members as they 
appear tn the fuel ertb array. Fire spread is from the left. 


least 0.6 cm. before coming through the back or bottom of the fuel member. Thus, the 
junction and the immediately adjacent leads lie in an isothermal surface minimizing 
heat transfer to or from the junction along the leads. The temperature field is 
captured by simultaneous measurements at the moment of ignition. (The scan time is 
extremely short compared to the time of the event.) Since only a limited number of 
thermocouples could be placed in the transverse fuel member, the highest density was 
concentrated in the region of the steepest temperature gradient, the upper corner 
nearest the advancing fire (fig. 4). Therefore, isotherms generated from the temper- 
ature field outside the upper corner are less reliable. 


Applying the same instrumentation procedure to the longitudinal fuel member 
would require an extensive array of thermocouples parallel to those located in the 
transverse members. The alternate method adopted was based on the average rate of 
spread of the thermal wave that travels along the longitudinal member. From the 
average rate of spread it is possible to convert temperature variations with time 
to temperature variations with distance. Thus, it is only necessary to determine the 
temperature variation with time at various prescribed vertical depths in the fuel 
member. The temperature distributions at each depth can be interrelated through the 
ignition time. Small holes were drilled to varying depths from the bottom (fig. 2). 
Chromel-alumel thermocouples were inserted the full depth of each hole and plugged 
with wooden dowels. There is no appreciable temperature gradient normal to the 
vertical plane determined by the thermocouples. The velocity of the thermal wave was 
obtained by superimposing the temperature histories of the two top thermocouples and 
noting their shift in time and separation. 
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Figure 4.--Cross-secttonal views of the locatton of thermocouples within the transverse 
and longttuditnal members of the untt cell. A ts transverse and B ts longttudinal. 
Thermocouple postttons are tndtcated by the crosses. 


The velocity of an advancing thermal wave in the longitudinal member, along with 


the rate of spread of the fire, is tabulated below. The packing ratio was held constant 
at 0.08. 


Stze Velocetty 
(Cm. ) Fire spread Thermal wave 
----------~ (Cm. /sec, )------------ 
0.6 0.33 20 .53* 
6 $5 5 Sie 
Wd 35 1.40 
lle 299 «95 


The discrepancy between the rate of spread of the fire and the rate of spread of 
the thermal wave within the unit cell is a result of the quasi-steady character of the 
fire. Viewed at close range, the fire appears to stop at each transverse member and 
then quickly travel down the longitudinal member to the next transverse member where 
it waits again. Viewed from a distance, the average rate of spread draws our attention 
away from the unsteady motion of the fire on the microscale. The thermal wave is 
probably more representative of the longitudinal fuel member than of the fuel array. 


2*Values obtained from adjacent longitudinal members in the same fire. 
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Results 


Isotherms obtained from longitudinal and transverse members of a unit cell within 
the fuel array are shown in figures 5 and 6, respectively. Two fuel sizes are repre- 
sented: 0.6 cm. and 1.3 cm. The transverse and longitudinal members are analyzed 
separately; no attempt is made to merge isotherms because no temperature data are 
taken in the region of the connection. Isotherms in this region are extrapolations. 
Consequently, the isotherms in the longitudinal members do not show masking by the 
transverse member near their intersection. 


Areal temperature distributions from figures 5 and 6 were used in equations (6a) 
and (6b) to obtain the nonuniform heat absorbed by the individual unit cell members up 
to ignition (table 1). No longitudinal member was instrumented during the first test 
at the 0.6-cm. sample size; however, two adjacent longitudinal members were instrumented 
in the second test at the 0.6-cm. sample size. 


Because final results are limited in figure 7, more than one function will fit the 
data. However, the exponential forn, 


e = exp (-4.53/c), 
is favored because: 
1.. The data lie close to the exponential curve; 


2. the curve point (e=1, 1/o = 0) acts as a constraint on the possible 
curves that can pass through the data points; and 


3. the region of the curve covering small particle sizes (near e=1) exhibits a 
linear relationship with 1/o, as suggested earlier by Thomas (1967) with 
A « 1l/o. 


Although data are confined to fuels small in cross section, extrapolation to 0.6 cm. is 


dependable. This judgment is based on reasonable results obtained by using the above 
exponential relationship in Rothermel's (1972) fire spread model. 
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Figure 5.--A cross-secttonal view of tsothermal surfaces withtn the transverse member 
of the unit cell at tgnttiton. Surfaces are parallel to the average combustton zone 
front. Each member ts tdenttfted bu test number for reference to table 1. 

B = 0.08. 
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Figure 6.--An averaged cross-sectional view of tsothermal surfaces within the longttu- 
dinal number of the untt cell at ignition. Surfaces are parallel to the average 
combustion zone front. Each member ts identified by test number for reference to 
table 1. 8B = 0.08. 
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Table 1.--Effective heating number. 


Test 3. article: “s ‘ : : G ; 
identification : thickness : Qn : QB : Qy*Qp : Tig : Qig : = 

----- Cal./g.- ---- nO eal. /g. 

1 0.6 108 1gqx 193 340 189 0.51 

2-R 6 94 aS 169 300 164 52 

2-L 0 94 FA 165 300 164 20) 

3 eS 75 19 94 380 216 s22 

4 LES) 72 20 93 330 183 BAS 


1 * Value is obtained by assuming that the ratio between the two values of Q 
for the 0.6-cm. specimens is equivalent to the ratio of the values of Q, where the 
known Qa is an average of the R and L (right and left of transverse member) values. 


2010 4 2 l 
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Figure 7.--A curve of €, 
the effective heating 
number, versus the re- 
ciprocal of o, the 
surface area-to- 

volume ratio. Stnce o 
= 4/t, the absctssa may 
be viewed as one-fourth 
the thickness of the 
particle. 
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Conclusions 


Nonuniform ignition over the unit cell weakens the analysis. This effect was 
minimized by choosing a small unit cell (8=0.08). The packing ratio here is represent- 
ative of the upper range of packing ratios occurring in natural fuels. 


The heating efficiency required to bring fuel ahead of a fire to ignition as it 
closes with the fire front has been investigated. Results show the dependence of ¢€, 
the fractional amount of the bulk density effectively brought to ignition, on the 
surface area-to-volume ratio of the fuel in the bed. At present, these data provide 
the only means for evaluating e. 


Structure and geometry have limited the number of sensing devices that can be 
installed in fuel members without greatly altering the measurement. Therefore, more 
thermocouples were placed in areas of greater sensitivity to the incoming flux. 
Combined results from the temperature profiles of the longitudinal and transverse fuel 
members provide a measure of the fractional amount of the bulk density that is involved 
in the ignition process. These results are reasonably sensitive to the surface 
area-to-volume ratio of the fuel member. 
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Nomenclature 
Definttton 
Fire spread rate 


Reaction intensity--the total intensity available 
from the combustion zone of a spreading fire 


Propagating intensity--that portion of the reaction 
intensity that drives the fire 


TR 


Effective heating number--the fraction of the bulk 
density that is effectively brought to ignition 
Ovendry bulk density 

Ovendry particle density 


Heat of preignition--the heat necessary to bring 
a unit ovendry mass uniformly to ignition 


Nonuniform heat absorbed per unit volume up to 
ignition 

Nonuniform heat absorbed per unit ovendry mass for 
the transverse member 


Nonuniform heat absorbed per unit ovendry mass for 
the longitudinal member 


Unit cell volume 

t?2--volume of transverse member of unit cell 

Square cross-sectional thickness of unit cell member 
Length of unit fuel cell member 

t2 

th 


Packing ratio--the ratio of occupied volume to 
the total bulk volume 


l--specific heat of water 

Co + C,T--specific heat of dry fuel 
0.266 

0.00116 

Temperature 

Ambient temperature 


Ignition temperature 


Fractional fuel moisture content relative to 
its ovendry mass 


Heat per unit dry mass absorbed by the transverse 
member of the unit fuel cell uniformly heated 
to ignition 


540--heat of vaporization of water 
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Units 


cm./sec. 
cal./cm.“-sec. 


cal. /cm.“-sec. 


g./em. 3 
g./cm. ° 


cal./g. 
cal./cm. 
ecalisice 


cal./s. 


cm. 3 


cm. 3 
cm. 
cm. 
cm. 


cm. 


cm. ?/em. 3 

Cal. /g7=2G. 

cal, /g.="C; 

Cali/ee=2Ge 
(cal./2.-"€.) 2G 

2G; 

vce 

nie 

g./g. 


cal 


Cal/ion, 
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Headquarters for the Intermountain Forest and 
Range Experiment Station are in Ogden, Utah. 
Field Research Work Units are maintained in: 


Boise, Idaho 

Bozeman, Montana (in cooperation with 
Montana State University) 

Logan, Utah (in cooperation with Utah 
State University) 

Missoula, Montana (in cooperation with 
University of Montana) 

Moscow, Idaho (in cooperation with the 
University of Idaho) 

Provo, Utah (in cooperation with Brigham 
Young University) 

Reno, Nevada (in cooperation with the 
University of Nevada) 
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